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Abstract 

The form factor for the anomalous process 77r* vrvr, F^'" {s,t,u), is calcu- 
lated as a phenomenological application of the QCD Dyson- Schwinger equa- 
tions. The chiral-limit value dictated by the electromagnetic, anomalous chiral 
Ward identity, ^^'^(O, 0, 0) = eNc/ (127r^/^), is reproduced, independent of the 
details of the modelling of the gluon and quark 2-point Schwinger functions. 
Using a parametrisation of the dressed u-d quark 2-point Schwinger function 
that provides a good description of pion observables, such as vr-vr scattering- 
lengths and the electromagnetic pion form factor, F^'^{s,t,u) is calculated 
on a kinematic range that proposed experiments plan to explore. Our result 
confirms the general trend of other calculations; i.e., a monotonic increase 
with s at fixed t and u, but is uniformly larger and exhibits a more rapid rise 
with s. 
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1. Introduction. Hadronic processes involving an odd number of pseudoscalar mesons 
are of particular interest because they are intimately connected to the anomaly structure of 
QCD. The decay vr*^ 'j'j is the primary example of such an anomalous process. That such 
processes occur in the chiral limit (m^ = 0) is a fundamental consequence of the quantisation 
of QCD; i.e., of the non-invariance of the QCD measure under chiral transformations even 
in the absence of current quark masses Jll. The 7r° 77 decay rate can be calculated from 
a quark triangle diagram and agreement with the observed rate requires that the number of 
colours, Nc, equal three. The transition form factor for the related process 7*7r° — >■ 7 can 
be measured experimentally |^ and has attracted keen theoretical interest |]^,^ because it 
involves only one hadronic bound state and provides a good test of QCD-based models and 
their interpolation between the soft and hard domains. 

Another anomalous form factor, accessible to experiment, is that which describes the 
transition 77r* — vrvr, denoted by F^'^{s,t,u). This provides additional constraints on QCD 
based-models because three hadronic bound states are involved. The form factor F^'"{s, t, u) 
has been measured at Serpukhov in the Primakov reaction tt^A ti^'ti^A'. [^] In this 
experiment the considerable uncertainty in both the kinematic range and result make it 
difficult to draw a conclusion regarding the accuracy of the theoretical prediction for the 
chiral limit value of ^^^^(0,0,0) 0. New experiments are planned at CEBAF: 0] 7P 
TT'^ir^n, s/m^ G [4, 15]; and at FermiLab p via the Primakov reaction using a 600 GeV pion 
beam, s/m^ G [4,6]. 

Herein we report a calculation of F^'^{s,t,u) on the kinematic range to be explored 
in the CEBAF experiment 0. We evaluate the amplitude indicated by the diagram in 
Fig. in which the quark 2-point Schwinger function, quark-photon vertex and quark- 
pion vertex (pion Bethe-Salpeter amplitude) are dressed quantities whose form follows from 
the extensive body of nonperturbative, semi-phenomenological Dyson- Schwinger equations 
studies in QCD. PJToO In this way our calculation provides for an extrapolation of the 
known large spacelike-g^ behaviour of these QCD Schwinger functions to the small spacelike- 

region, where they are unknown and confinement effects are manifest. This facilitates 
an exploration of the relationship of physical observables to the nonperturbative, infrared 
behaviour of these Schwinger functions. 



This calculation employs the model forms introduced in Ref. |T^. The quark 2-point 
Schwinger function has no Lehmann representation and hence may be interpreted as de- 
scribing a confined particle since this feature is sufficient to ensure the absence of quark 
production thresholds in S'-matrix elements describing colour- singlet to singlet transitions. 
The quark-photon vertex, which describes the coupling of a photon to a dressed quark, 
follows from extensive QED studies |]TT],|T^ and satisfies the Ward-Takahashi identity. This 
necessarily entails that the amplitude is current conserving. In the chiral limit the quark- 
pion vertex is completely determined by the quark 2-point Schwinger function [|l3l, which 
is the manifestation of Goldstone's theorem in this approach. The extension to finite quark 
mass requires a minimal modification and preserves Dashen's relation [|14| . 

The quark 2-point Schwinger function is parametrised such that it provides an good 
description of vr-vr scattering at, and to approximately 600 MeV above, threshold and the 
electromagnetic pion form factor at spacelike momentum transfer. No adjustment of the 
parameters is made in calculating F^'^{s, t, u). The 77r* — vrvr process probes the pion bound 
state amplitude well outside the domain of the complex plane on which it has been fitted; 
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i.e., well into the timelike region, and hence provides a new test of the model. Therefore, 
interpreted within our framework, the experimental determination of F^'"{s,t,u) provides 
important new information about the structure of the nonperturbative pion bound state 
amplitude. 

2. The amplitude for 777* — > nn. In Euclidean space, with metric 
6fj_u = diag(l, 1, 1, 1) and 7^ = 7^^, the amplitude for the 77r*7r7r vertex (see Fig. |lD is 

2eA^c/ (0trz, [r^(A:___,A;+++)5(A;+++)75r.(W)5(A:-++) 
75r.(A:-o+)5(fc-_+)75r.(A;__o)5(fc— )] • (1) 



In this expression the colour and flavour traces have been evaluated, leaving only the Dirac 
trace, and we have employed the definition = k + + |p2 + fPs- Note that due 
to momentum conservation the photon momentum is equal to the sum of pion momenta: 
q = Pi + P2 + Ps- In the following we adopt the convention that the pions labelled 1 and 
2 are on-shell whereas pion 3 is off-shell. The dressed quark-photon vertex is denoted by 
1^/^(^15^2)5 the pion Bethe-Salpeter amplitude by 75r^(/c) and the dressed quark 2-point 
Schwinger function by S{k). Equation (0) can be derived as an application of the formalism 
described in Ref . . Evaluated with a dressed quark-photon vertex that satisfies the Ward- 
Takahashi identity, which is a minimal requirement, this expression is current conserving. 
The quark 2-point Schwinger function can be written 

S{p) = -ij-pav{p^) + crs{p^) (2) 

= ^ , (3) 

i'j ■ p A{p'^) + m + B{p^) ' 

where m = m„ = is the current quark mass, and can be obtained by solving the quark 
Dyson-Schwinger equation (DSE) |p. The many studies of this equation ensure that the 
qualitative features of the functions as and ay are well known for real, spacelike-p^ . 

In Ref. |]TU[, in order to avoid the need for a numerical solution of the quark DSE, the 
following approximating algebraic forms are used: 

-as{x) = Cmc-'- + L + b2^-^] + (1 - e-(-+-^)) (4) 

OiX 03X Y Ax J X + m"^ ^ ' 

l[x -|- m^Y 

with k"^ = 2Dx, (fsix) = \/2D asik"^), o^vix) = 2Dav{k'^) andm = m/\/2D. The quantity 
is a mass-scale related to the infrared behavior of the gluon 2-point Schwinger 
function 0. 

At large spacelike-p^ the model forms behave as 



cTsip^) ~ - - — + riT- + ■ ■ ■ ^^(P ) ~ — + • • • ' 6 

p^ 01O3 p^ p^ p^ 

whereas in QCD one has 
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<^s{p') ~ — - r / + • • • (^) 

p2[lln(p7A2jCD)] ^ p^[ln(pVA2^cD)] 

with m and (gg) the renormahsation point invariant current mass and condensate, respec- 
tively, and d = 12/ [33 — 2Nf]; Nf is the number of quark flavours. Comparing these 
two equations one sees that, setting d = 1 and neglecting ln[p^] terms, the model defined by 
Eqs. (ID and (|^) properly represents the ultraviolet behaviour of the quark 2-point Schwinger 
function; incorporating asymptotic freedom and dynamical chiral symmetry breaking, with 

Another feature of this model is that both as and cry are entire functions in the complex- 

plane but for an essential singularity. As a consequence the quark 2-point Schwinger 
function does not have a Lehmann representation and can be interpreted as describing a 
confined particle. This is because, when used in Eq. (|I|), for example, this property ensures 
the absence of free-quark production thresholds, under the reasonable assumptions that 
Tt^Ip^) is regular in the domain of integration. It follows from this that Eq. (|l]) is free of 
endpoint and pinch singularities. 

The expressions in Eqs. (H) and (|^) provide a six-parameter model of the dressed quark 
2-point Schwinger function in QCD: Cm, fh.bQ, . . . , 63. (A = 10""^ is introduced simply to de- 
couple 62 from the quark condensate.) These parameters can easily be fitted to experimental 
observables, as we discuss below. 

The Bethe-Salpeter amplitude, F^r in Eq. (|l|), is the solution of the homogeneous Bethe- 
Salpeter equation (BSE). Many studies of this BSE suggest strongly that the amplitude is 
dominantly pseudoscalar (as already indicated by the notation used in Eq. ([1|)). Further- 
more, in the chiral limit the pseudoscalar BSE and quark DSE are identical [|13| and one has 
a massless excitation in the pseudoscalar channel with 

F^(p;p2 = o) = -^i?^=o(p') , (9) 

where Bm=o{p^) is given in Eq. (|^) with m = 0. This is the realisation of Goldstone's theorem 
in the DSE framework; i.e., in the chiral limit Eqs. (^) and (^ completely determine F^. 
Herein we employ the approximation 

T^{p-P' = -ml)^^Bm=,{p') , (10) 

which, for small current quark masses, is a good approximation both pointwise and in terms 
of the values obtained for physical observables. fl^ 

The quark-photon vertex, F^(pi,p2), satisfies a DSE that describes both strong and 
electromagnetic dressing of the vertex. Solving this equation is a difficult problem that 
has only recently begun to be addressed |T^. However, much progress has been made in 



constraining the form of F^(pi,p2) and developing a realistic Ansatz |11,12|. It is obvious 
that the bare vertex, F^(pi,p2) = is inadequate when the fermion 2-point Schwinger 
function has momentum dependent dressing because it violates the Ward-Takahashi identity. 



In Ref. IjT^ the following form was proposed 
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2 

{P + k)f. 



7m 



p2 _ J^2 



[7 ■ p + 7 ■ /c] 



— I 



This Ansatz is completely determined by the dressed quark 2-point Schwinger function and 
has the features that it: 1) satisfies the Ward-Takahashi identity, thereby ensuring current 
conservation at the microscopic level; 2) is free of kinematic singularities - i.e., it has a well 
defined limit as p"^ k"^; 3) transforms correctly under C, P, T and Lorentz transforma- 
tions; and 4) reduces to the bare vertex in the manner prescribed by perturbation theory. 
The fact that it is nevertheless relatively simple makes it an ideal form to be employed in 
phenomenological studies. The studies of Ref. [|lOl illustrate its phenomenological efficacy. 
In these studies, for example, it is shown that the chiral limit for the anomalous decay 
^ 77 is reproduced exactly using 



TT 



(12) 



3. Chiral Limit: 777* — > tttt. At the soft point in the chiral limit {s = t = u = 0) the 
transition form factor of Eq. (P is 



F3-(0, 0,0) = ^ ds s T^{s f ay [A oy as + \sA a'y as - \sA ay a's 

+ \sA! ay as - \sB' a\^ 



(13) 



Defining C(s) = Bi^sf' l\s A{sY\ = crs(s)^/[s ^^/(s)^] one obtains a dramatic simplification 
and Eq. (0) becomes 



^3^(0,0,0) 



eiV, /-^^ r^(s)3 C\s)C{s) 



27r2 Jo B{sf \\^C{s)\ 



(14) 



Recalling Eq. ( 
it follows that 



), which is the manifestation of Goldstone's theorem in the DSE approach. 



F=^^(0,0,0) 



dC- 



C 



1 + cy 



(15) 



since C{s) is a monotonic function for s > with C(s = 0) = cxo and C{s = 00) = 0. 
Hence, the chiral limit value is reproduced independent of the details of the quark 2-point 
Schwinger function, S{p). 

In order to obtain the result in Eq. (|15|) it is essential that, in addition to Eq. (^, the 
photon-quark vertex satisfy the Ward identity. This is not surprising. However, the fact 
that one must dress all of the elements in the calculation consistently is often overlooked. 
The subtle cancellations that are required to obtain this result also make it clear that it 
cannot be obtained in model calculations where an arbitrary cutoff function (or "form- 
factor" ) is introduced into each integral. The fact that the pion Bethe-Salpeter amplitude is 
proportional to the scalar part of the quark self energy in the chiral limit, Eq. (^, is crucial. 
These features are also to be seen in the calculation of the Wess-Zumino five pseudoscalar 



term and the vr 77 vertex |jTO[, in which again the values expected from considerations of 
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anomalous current conservation are obtained independent of the details of the quark 2-point 
Schwinger function, S{p). 

4. Numerical Results. The chiral limit value of the '-j-k-k-k amplitude, F^'^(0,0,0), 
provides a useful normalisation. We therefore define the function 

F'-{s,t,u) = F-'-ip,,p2,ps)/F'-i0,0,0) = Un'f^F'''{p,,p2,P3) . (16) 

We employ the following definition of the Mandelstam variables: 

5 = -{pi + P2f = mis, t = -pi = mli, u = -{pi + p^^f = mlu, (17) 

which ensures that even though we work in Euclidean metric these variables have their 
conventional interpretation. We note that the quantity t — provides a measure of the 
amount by which the (third) pion is off-shell. 

In the experiment proposed at CEBAF the photon energy in the proton rest frame is 
between 1 and 2 GeV. This suggests the following range of Mandelstam variables: 4 < s < 
16, —9 <i< —1, —16 < M < 5. For the purpose of the numerical calculation we fix t = — 1; 
i.e., we choose pion 3 to be as close as (experimentally) possible to its mass shell, t = 1. 
Within the range of s considered, the requirement of a fixed photon energy of 1 GeV in the 
proton rest frame entails 

M = (-1.5 -4.28 x + 0.012 a;2) with x = (s/4 - 1). (18) 

We have demonstrated above that in the chiral limit F'^'^(0,0,0) = eNc / {12tt'^ f^) in- 
dependent of the model parameters. The evolution of F^'^{s,t,u) with s does depend on 
the model parameters, even in the chiral limit. As we have described above, Eqs. (^, 
(I), (H) and ( |T2|) provide a six-parameter model [Cm, rn, 6o • • • ^sjof the nonperturbative. 



dressed-quark substructure of the pion based on DSE studies. These parameters are fixed 
by requiring that the model reproduce, as well as possible, the following pion observ- 
ables U/{qq)l^^^y2 = 0.42 ± 0.02, f^r^ = 0.31 ± 0.01, ml/{qqf/^^^, = 0.40 ± 0.03; the di- 
mensionless vr-vr scattering lengths (discussed in Refs. pO| , ^^ with current experimen- 
tal values presented in Table |) and the pion electromagnetic form factor on spacelike- 
e [0,4] GeV2 0. 

The values of observables are given by simple integral expressions involving the quark 



2-point Schwinger function and pion Bethe-Salpeter amplitude; for example, [!(: 



= (19) 

^ dp" p" Bo{pY ((t^ - 2 [asa's + saya'y] - s 



87r 



2 ^2 

mt JZ 



27r2 



. 2 2 Boip"^ 
5(^ 



{B{p')a%s)-Boip')crsip' 



(20) 



which follows from the pion Bethe-Salpeter equation and is consistent with Dashen's relation 
I^ ; and {qq)^2 = In fi'^/AqQ^) {qq) with {qq) given in Eq. and Aqcd = 0.2 GeV. In 
these equations the subscript or superscript "0" indicates that the labelled function has 
been evaluated with zero quark current mass. 
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Following this procedure one obtains 

Co = 0.121, m = 0.00897 , Cm = 0, 

bo = 0.131 , bi = 2.90 , 62 = 0.603 , 63 = 0.185 . (21) 

The mass scale is set by requiring that take its experimental value of 92.4 MeV, which 
yields D = 0.160 GeV^. The calculated values of observables are presented, for completeness, 
in Table |. The form factor is discussed in Refs. [10 

The amplitude in Eq. (P is given by a four- dimensional integral that we evaluate by 
straightforward Gaussian quadrature. The calculation is simplified by working in the rest- 
frame of one of the outgoing pions. All results reported here are tested against variation of 
numerical parameters, such as grid densities, etc. 

The quark 2-point Schwinger function specified by Eqs. (H), (H) and (|) is an entire 
function. The quantity x[(Ty(x)]^ -|- [^■^(x)]^ has a pair of complex conjugate zeros at 
x = —0.102 ± i 0.715, which corresponds to {p^ = x/{2D) \ ~ 12 m^. This affects our calcula- 
tion through the approximation of Eq. (^) because -Bo(p^) has poles at these points. These 
poles are integrable singularities and do not lead to an imaginary part in the amplitude. 
However, in our calculation of F^'"{s,t,u) we treat them as a spurious artifact of the ap- 
proximation of Eq. ( p!OD and consider our results to be unreliable when these poles enter the 
integration region. This occurs only for s > 12. The approximation of Eq. ([lO|) has hitherto 
only been tested for real, spacelike-p^ and the present application involves a considerable 
extrapolation into the complex-p^ plane as s is increased. (In the chiral limit, = 0, the 
integration does not explore the complex-p^ plane and these poles do not enter the integra- 
tion region.) The importance of this observation is that processes such as 77?* — vrvr provide 
a means of exploring the structure of bound state amplitudes in the complex-p^ plane. 

Our numerical results are shown in Fig. ^ The solid line is our calculated result for 
the 77r*7r7r amplitude F^'^{s,t,u), Eq. (|TBp, as function of the Mandelstam variable s for 
different values of u as given by Eq. ([T8| ) at i = —1. An excellent fit to this curve is given 
by 

F^''(s,t,u) = 1.044 + 0.096 a; + 0.006 with a; = (s/4-l). (22) 

Fixing the Mandelstam variables such that the energy of the incident photon is 2 GeV in 
the proton rest-frame changes our curve by an amount that is not visible on the scale of this 
plot. Importantly, the result is insensitive to the details of the parametrisation of the quark 
2-point Schwinger function. A quantitatively similar curve is obtained using earlier sets of 
the parameters in Eq. (pi]); i.e., our result is not sensitive to the details of the model. The 
form factor does depend on the pion mass. The short-dashed line is the result we obtain 
with = m^^P/2. 

A comparison of our result with that obtained in other models reveals that all results 
are broadly consistent. Our prediction for this form factor is, however, uniformly larger 
and displays a more rapid increase with s. For example, using vector meson dominance one 
obtains, for real photons p6[1 , 

F'-is, t,u) = l + C,e^' + + 
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with f = s — t + u — 2ml, = 2fi'p7r7rfi'p7r7/["^pF^''(0, 0, 0)] ^ 0.434 and 5 an unknown phase. 
This expression, with 5 = 0, appears as the long-dashed hne in Fig. This curve is smaller 
in magnitude than our result for all s and rises more slowly with s. Another vector meson 
dominance model estimate of this process appears as the dash-dot line in Fig. ||. This 
result is obtained from 



S, t, u) 



+ 



1 



V 



+ 



— u ^ 



(24) 



and is even smaller in magnitude and has a weaker s dependence than the other vector 
meson dominance estimate. 

The weakest s dependence is obtained using a model based on chiral expansion techniques 



and employing a vector meson saturation Ansatz |^ . This result appears as the dotted line 
in Fig. and is obtained from the expression 



1 ~l~ C'pion— loops 



S + t' + U 



(25) 



where the non-divergent part of the coefficient Cpion-ioops is 



pion— loops 



with /(C < 0) defined by 





'ml' 




^log 








ml_ 





(26) 



/(C) = (l-C/4m2).yi-4m2/Clo, 




4ml/ C + 1 



- 2 . 



4m2/C-l, 



(27) 



Analytic continuation is used to define f{( > 0). Note that f{() has an imaginary part for 
( > 4ml; ^ domain explored by the existing and proposed experiments. The 

imaginary part is due to the pion loop in the s-channel. 

We can compare our result with the one existing data point |^. Its statistical and 
systematic error as well as the uncertainty in s is also displayed in Fig. The fact that this 
data point is well above the chiral limit prediction has caused some concern However, 
given the experimental errors and the prediction of our model this data point does not 
appear untenable. 

5. Summary and Conclusions. Herein we have reported a calculation of the form 



factor for the anomalous process 77?* 



vrvr, 



F^'"{s,t,u), in a phenomeno logical approach 



based on the QCD Dyson-Schwinger equations (DSEs). In our approach the chiral limit 
value of ^^'^(O, 0, 0) is reproduced independent of the details of the quark 2-point Schwinger 
function. Using for the w-d-quark 2-point Schwinger function a parametrisation fitted to 
low-energy pion data, F^'^{s,t,u) was calculated for a range of {s,t,u) specified in such a 
way as to cover the kinematic region to be explored in a proposed experiment [0 . The small 
photon virtuality in another proposed experiment ||^ should also make it possible for our 
prediction to be compared, without adjustment, with the results of that experiment. 

We compared our result with that obtained in other models when applied in the same 
(s, t, u) range as we have considered, which is the range appropriate for the experiment 
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proposed at CEBAF [0]. Our model yields a form factor, F^'^{s,t,u), that is uniformly 
larger and has a more rapid increase with s than any of the other models considered. The 
result obtained in a model based on chiral expansion techniques, augmented by vector meson 
saturation assumptions, leads to the weakest s dependence. Although the models are broadly 
consistent, the variation between the results is such that the proposed experiments should 
be able to distinguish between them. 

Along with the study in Ref. this phenomenological application of the QCD DSEs 
is one of the first to explore the model quark 2-point Schwinger function and pion Bethe- 
Salpeter amplitude in the timelike region, which is not accessible in perturbation theory. 
This region is important in the study of, for example, vector meson Bethe-Salpeter and 
baryon Fadde'ev amplitudes. Experimental data on F^'^{s,t,u) can therefore be used to 
place additional constraints on the analytic structure of the QCD Schwinger functions. 
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TABLES 





Calculated 


Experiment 




0.0924 GeV 


0.0924 ± 0.001 




0.1385 


0.1385 


1 rpv2 


0.0051 


0.0075 


-(aa)^ 


0.221 


0.220 


' TT 


55 fm 


n _|_ n nna 

\J •\J\JtJ —1— \J •\J\J\J 




0.505 (dimensionless) 


0.504 ± 0.019 




1.04 


1 (Anomaly) 




0.17 


0.21 ± 0.01 


al 


-0.048 


-0.040 ± 0.003 


al 


0.030 


0.038 ± 0.003 


aO 


0.0015 


0.0017 ± 0.0003 




-0.00021 





TABLE L Pion observables calculated using the parameter values in Eq. (21). The "experimen- 
tal' values listed for m^™ and {qq) are an indication of other contemporary theoretical estimates. 
Experimental values not discussed in the text are taken from Ref. ||2^. The difference between 
the calculated and experimental values of is a measure of the importance of final-state tt-tt 
interactions and photon- p-meson mixing |24|; that between the calculated and experimental values 
of the pion scattering lengths is a measure of the importance of tt-tt final-state interactions in this 
case p0|. 
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FIGURES 



TT 



Pa 



7 




TT 



TT 

FIG. 1. This figure is a pictorial representation of the amplitude identified with the 77r*7r7r 
vertex. The straight broken external lines represent the outgoing vr, the circles at the vr legs 
represent the (7r|gg) Bethe-Salpeter amplitudes, F^, the wiggly line represents the photon, 7, the 
circle at the 7 leg represents the regular part of the dressed quark-photon vertex, F^, and the full 
internal lines represent the dressed quark 2-point Schwinger function, S. 
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FIG. 2. The normalized 77r*7r7r amplitude F^'^{s, t, u), Eq. (^), as function of the Mandelstam 
variable s for different values of u as given by Eq. (18) at t = —1. Our result: solid line; our 
result with ttt-tt mTr/2: short-dash line; vector meson dominance: long-dash and dash-dot lines, 
see Eqs. (|23| ) and (p^; chiral expansion plus vector meson saturation Ansatz: dotted line, see 
Eq. (p5[). The data point is taken from Ref. B. 
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